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I. INTRODUCTION 


Ultra-High Performance Fiber 
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Reinforced Concrete 
(UHPFRC) is an innovative material that can reach 
average compressive strengths at 28 days that surpass 
150MPa (22ksi), with tensile strengths of 7MPa (1ksi), and 
10MPa (1.5ksi) in bending. To obtain a mix with ultra- 


Abstract—Ultra-high performance fiber reinforced concrete (UHPFRC) 
is an advanced composite material characterized by compressive and 
tensile strengths above 150MPa and 7MPa, respectively. Initially, an 
experimental procedure was used to characterize the tensile performance 
through bending tests, using beams with 1% and 2% content by volume of 
steel fibers. Three-point bending load arrangement notched prisms were 
used to determine the contribution of the fibers to reinforcing a cracked 
section. With that, the (F vs. œ) experimental curves were graphed, and 
from there, the analytic tensile curves (o vs. œ) was obtained point by 
point by application of the inverse analysis procedure proposed by the 
AFGC. With the analytic curves, the fracture energy was calculated, 
following a procedure proposed by RILEM. Subsequently, the crack width 
was transformed into strain using a relationship that involves the 
characteristic length. The resulting analytical behavior law was used to 
carry out computational modeling applying the finite element method. 
Both the finite element method and the fracture energy were used to 
validate the procedures, comparing experimental and numerical results. 
Models and experiments showed good agreement and finally was 
determined the constitutive law for the UHPFRC in tension. It can be 
concluded from this study, therefore, that the post-cracking tensile 
behaviour of UHPFRC can be appropriately evaluated and validated 
through the applied analysis procedure in this research. 


high-strength, Camacho E. [1] observes that the water 
amount not chemically combined with the cement in the 
hydration process to be the less as possible, minimizing 
porosity and its connectivities, and increasing strength and 
durability. Schmidt and Fehling [2], additionally, have 
indicated that the particle packing should be optimized by 
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using large amounts of superplasticizers, adjusting the mix 
workability with the presence of fibers. 


Four principles that must be met to achieve ultra-high 
strength and durability in concrete: (i) a very low 
water/cement ratio, which in our case was 0.19, resulting 
in a very dense and strong structure, minimizing pore 
capillarity and preventing the transport of toxic gases and 
liquids into and through the concrete; (ii) high particle 
packing, this requirement was not fulfilled, in our case a 
simple grinding process was carried out.; (iii) the use of a 
large amount of superplasticizer, to adjust workability; (iv) 
the use of fibers to increase tensile strength, flexural 
strength, and shear strength and to make the concrete 
sufficiently ductile. By keeping such general design rules, 
it is possible to define UHPFRC mixes for the use in 
beams, which should continue to have a great deal of 
bending strength even after cracking. During the post- 
cracking behavior, the fibers, subjected to tensile, have a 
fundamental role, since, when appropriately oriented, they 
tend to prevent a fragile rupture due to their bridging 
action that sews both sides of cracks together. 


Naaman and Reinhardt [3] came up with a 
classification for Reinforced Concrete Fiber (FRC) that 
can be applied to UHPFRC. They classify FRC 
accordingly to the inelastic behavior: (i) tensile strain 
hardening, where the maximum internal force in the 
cracked zone is larger than the limiting elastic internal 
force, and (ii) tensile strain softening, where the limiting 
elastic force is larger than the maximum internal force. 
The recommendations made by the Association Frangaise 
de Génie Civil AFGC [4] indicate that when UHPFRC is 
subjected to the tensile, it can present both behaviors just 
mentioned, as well as define an intermediate behavior, 
named as low strain hardening. In this work, UHPFRC 
beams with steel fiber content by volume of 1% and 2% 
were subjected to three-point bending tests in the lab. 
Their responses, in terms of load vs. deflection (F vs. ò), 
were recorded and showed herein graphically. 


From there, the analytic tensile curves (o vs. œ) was 
obtained point by point by application of the inverse 
analysis procedure proposed by the AFGC [4]. With the 
analytic curves, the fracture energy was calculated, 
following a procedure proposed by the International Union 
of Laboratories and Experts in Construction Materials, 
Systems and Structures, RILEM TC50 [5]. The crack 
width was transformed into strain using a relationship that 
involves the characteristic length. The resulting analytical 
behaviour law was used to carry out computational 
modeling applying the finite element method. The program 
ANSYS [6] was used to carry out computational modeling 
and obtain the analytical load vs. deflection curves. This 
program requires, as input data, the constitutive behavior 
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of the material in compression and in tensile. The former 
can be specified with values directly obtained from the lab 
experiments, while the latter can be set with the Inverse 
Analysis. Both the finite element method and the fracture 
energy were used to validate the procedures, comparing 
experimental and numerical results, see Fig. 1. 
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Fig. 1: Research scheme 


Il. EXPERIMENTAL PROGRAM 


A practical strategy widely used in experimental programs 
to analyses the results of concrete resistance tests is the 
factorial arrangement, in which different treatments that 
are to be compared are defined. In the design of 
treatments, controllable factors, their levels and the 
combination between them are selected. The experimental 
design indicates the way in which the treatments are 
randomized and the way to control their natural variability. 


The statistical tools indicated previously were used to 
define the UHPFRC mixture design used in this study, 
which was part of a series of studies carried out as support 
to an invention patent. It was deposited with the National 
Institute of Industrial Property (INPI), see Rojas et al. [7]. 
Also, those developments can be consulted in more detail 
in Rojas R. [8], it explains the extensive experimental 
work carried out, the end result of which is the design of 
the mixture indicated in TABLE 1, allowing the 
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production of UHPFRC with a compressive strength 
greater than 150 MPa. 


Table 1: UHPFRC mix design. 


Material kg/m? 
Cement 955 
GGBS 263 
Silica Fume 119 
Quartz powder 119 
Fine sand 788 


Superplasticizer 40 
Water 185 


2.1 Materials 
The agglomerating materials used in the mixture are made 
up of: 
— National cement type Portland CP V ARI with high 
initial resistance. 


-Ground Granulated Blast Furnace Slag (GGBS) 
donated by the company ArcelorMittal Tubarão in the 
Brazilian state of Espirito Santo. In TABLE 2 we 
observe the chemical composition of the GGBS used in 
this research, including the ranges recommended in 
ACI 233R [9]. 

— Commercial silica fume (SF). 

— Commercial quartz powder. 

It has a single aggregate consisting of silica sand with a 
maximum grain size of 0.30 mm. A solution of 
polycarboxylate in an aqueous medium (Visco-Crete 3535) 
supplied by SIKA was used as a super-plasticizer additive, 
which adjusts the workability of the concrete and is mixed 
with normal water to be placed in the mix. 


Table 2: Chemical composition of GGBS 


Main chemical Percent by Range ACI 
constituents mass [25] 
CaO 44.50% 32-45% 
SiO2 30.22% 32-42% 
A1203 7.92% 7-16% 
Fe203 7.45% 0.1-1.5% 
MnO 1.10% 0.2-1.0% 
MgO 1.08% 5-15% 


The fiber used is of the steel Dramix type, 13 mm long 
and 0.2 mm in diameter, in a volume equal to 1%. TABLE 
1 shows the proportions of the mixture, in which 26% of 
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the cement is replaced by sustainable materials (GGBS and 
SF) and 8% is replaced by quartz powder. The 
water/cement ratio is 0.19 and the water/binder ratio is 
0.13. 


2.2 Manufacturing the Mixture 


A sustainable UHPFRC is produced in this research. It has 
a simple manufacturing process, without the need for 
elaborated and delayed grinding processes for the 
packaging of particles. Two types of industrial waste are 
included in the mixture, silica fume and mainly GGBS, the 
latter with a specific granulometric distribution indicated 
in the invention patent. 


The materials are weighed and placed in a mixer in the 
following order: silica fume, cement, and blast furnace 
slag and silica sand. The dry materials are mixed for about 
5 minutes before the superplasticizer previously mixed 
with the water is added to the mixture. Wet materials are 
mixed for about 10 minutes. Initially, a dry mix is 
observed until small spheres of material are formed; about 
1 mm in diameter, these spheres get mixed together and 
progressively increase in diameter until they become a wet 
concrete paste. 


It is observed how the material separates from the 
bottom of the mixer, acquiring the shape and consistency 
of a dense plastic mass, see Fig. 2. In this state, the mixture 
for the UHPFRC is considered ready and it is in this 
moment that the steel fibers are placed, mixing for 
approximately 2 minutes. After fabrication, the mixture is 
cast into the respective moulds, to be compacted on a 
vibrating table for 1 minute. The specimens are stored and 
covered with a plastic layer for 48 hours, after which they 
are placed in a thermal bathroom for 24 hours at a 
temperature of 60 °C and then at 90 °C for another 24 
hours. They are then stored in a humid room at 23 + 3 °C 
until the day of the test, avoiding in all cases thermal shock 


on the specimens. 


Fig. 2: Mixture consistency. 


2.3 Post-peak uniaxial compression test 


From an experimental point of view, the compilation of 
consistent and accurate stress vs. strain data (o-s) is 
difficult. During the execution of the compression test, 
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when the first crack forms, the lateral strain exceeds its 
tensile capacity and the UHPFRC specimens (with fibers) 
behave elastically up to approximately 80 to 90% of their 
compressive strength. 


After reaching the maximum resistance (fc), a 
progressive strain softening takes place in which the 
presence of fibers regulates the softening stage in a similar 
way as it happens in tensile, to later produce the ductile 
compression failure. Hassan et al. [10] found a post-peak 
measurement method, which consists of placing the 
circular rings with the LVDTs in the specimen only to 
measure the elastic state of the test. Additionally, two 
LVDTs are placed parallel to the specimen to measure the 
movement of the test machine head, allowing the 
recording of the post-peak stage. That method was used in 
this research to recording the post-peak behaviour of the 
UHPFRC subjected to uniaxial compression. 


The uniaxial compression test was performed on 
specimens manufactured using steel moulds of 50mm 
diameter by 100mm height, containing a 1% fiber volume 
and following the criteria specified in the ABNT 
NBR7215 standard [11]. Twenty specimens, with 28 days 
of cure, were tested, applying monotonic displacement 
loading, using a 2000KN hydraulic machine at a rate of 0.5 
mm/min. Previously, the superior and inferior face of each 
cylinder was levelled mechanically using a rectifier and its 
height is measured to verify the necessity of applying 
some correction factor in the resistance according to 
ABNT NBR5739 [12]. 


The values of the load vs. vertical displacement of each 
specimen are recorded. In the linear elastic part, the value 
of the strain is calculating by dividing the average 
displacements of the LVDTs by the initial length of 
measurement maintained by the circular rings. Later, with 
the appearance of the first crack, a multiple cracking phase 
occurs, in which the strain is obtained by dividing the 
average displacement of the external LVDTs (those that 
measure the displacement of the machine head) by the 
total height of the specimen. The stress in this stage was 
obtained by dividing the machine load by the cross- 
sectional area of the cylinder. 


The characteristic compressive strength (fck) of the 
UHPFRC was calculated using the AFGC [4] 
recommendations, and the following considerations were 
taken into account: 


— Apply the displacement control load. 
—The specimen must exhibit a conical failure pattern. 


—The average strength must be calculated on at least 
three specimens. 


—The characteristic compressive strength value must be 
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calculated by subtracting the Student's coefficient 
multiplied by the standard deviation from the average 
strength value. 


2.4 Modulus of elasticity 


The modulus of elasticity was calculated by measuring 
directly on the linear upward branch of the UHPFRC 
constituent curve, recorded for each of the uniaxial 
compression tests performed on cylindrical specimens. 


A linear approximation is used with best fit o-z results 
between 0 and 80 % of the peak compression strength. The 
value of Ecm is then defined as the average modulus of 
elasticity of the UHPFRC or the average secant modulus 
of elasticity, calculated as the average of the twenty 
individual values obtained graphically. 


2.5 Three points bending test 


Three-point bending load arrangement notched prisms 
were used to determine the contribution of the fibers to 
reinforcing a cracked section. With that, the (F vs. œ) 
experimental curves were graphed. Ten beams (four with 
1% of fiber content and six with 2%), were manufactured 
with the mix presented in TABLE 1 and with the 
dimensions of 10x10x40cm. The lab tests were carried out 
in a hydraulic universal testing machine with a capacity of 
2000KN (450kip), after 28 days of curing and by applying 
displacements at a speed of 0.5mm/min (0.02in./min). All 
of them had a notch of 30mm in depth by 4mm in width at 
the bottom centre of their span length made with a circular 
saw. 


A horizontal LVDT type sensor was placed to measure 
the opening of the notch (w) and two vertical LVDTs, 
placed on each side of the beams, were used to measure 
their central deflection (6), as illustrated in Fig. 3. 


Fig. 3: LVDTs to measure œ and ô. 


The sensors are attached by means of tabs glued, a fast- 
setting glue is used. The distance between tabs must be the 
same from one test to the next so that the initial 
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measurements can be corrected by subtracting the elastic 
strain. 


The distance between tabs should be less than 4 or 5 
cm and the stroke of the sensors must be at least 2 mm. 


The test variability was controlled by using materials 
from the same batch and the same equipment to 
manufacture and test the specimens. The data results were 
digitally recorded for each test and graphically analyzed in 
load vs. deflection (F vs. 5) curves, as well as in load vs. 
crack width (F vs. œ) curves for each of the beams. 


NI. EXPERIMENTAL PROGRAM 


The relationship between forces and strains can be directly 
determined when: the internal force in a structural element 
is uniaxial; the cross-section is known, and it is possible to 
directly measure the deformation on the element under the 
action of a load. 


Using experimental data from uniaxial compression 
tests, o vs. € curves were obtained for the UHPFRC, which 
were then used as a part to input data for the computational 
modeling. 


When the forces were in bending, the determination of 
the nonlinear strains was not direct and alternate analytical 
procedures had to be used in the calculations. The 
procedure to determine the constitutive law for the 
UHPFRC in tension, including the post-cracking response, 
followed the methodology by AFGC [4]. 


The tensile curve was obtained point by point by 
application of the Inverse Analysis, i.e., obtaining the o vs. 
e analytic curve from the F vs. œ experimental curve. 


Both curves o vs. € in compression and in tension were 
introduced as input data for the computational modeling 
and then the F vs. 6 analytical curve was obtained. 
Therefore, a graphical comparison between the 
experimental and the analytical behaviors for each of the 
specimens tested were carried out. 


3.1 Inverse analysis 


The process starts with the definition of a new coordinate 
system at the point where the first crack occurs. The notch 
opening value at that point is turned into the new origin, 
with the first point coinciding with the elastic limit. 


The equilibrium is easily solved to find the internal 
force. From the first point (step i), the next points are 
calculated (steps i+1) by solving the equilibrium of the 
cracked section. 


A complex nonlinear equation system is generated at 
each step and, therefore, the free software Maxima [13] 
was used as a mathematical tool to solve the equations. 
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After solving the equation system, the force at the point 
is calculated, i.e., in this case, the cohesive force. The 
process is repeated at each i+1 point until the curve of 
cohesive force versus notch opening is built (actually, the 
oc vs. œ curve). Then, the oc vs. œ curve is transformed 
into a o vs. € curve, which, according to AFGC [4], can be 
used to define a relation between œ and e mainly based in a 
determination of the characteristic length (Ic) see equation 
1: 

fct,el w 


where: 

fct,el is the tensile strength of the concrete matrix; 

E is the modulus of elasticity of the concrete matrix; 
Ic is the characteristic length; and 

@ is the notch opening. 


The characteristic length is measured at the location 
where cracking occurs and in the same direction of the 
bottom notch opening of the beam. In the case beams are 
subjected to three-point bending, the AFGC [4] defines the 
Ic value as a function of the type of experimental behavior 
that is presented, i.e., the value depends upon the behavior 
as either of the strain softening or strain hardening types. 


If the beam presents a strain softening type of behavior, 
the characteristic length is calculated with equation 2, 
while if presenting a strain hardening behavior, equation 3 
is therefore used. 


= 2009 (2) 
E. G 
Ic= Er (3) 


where: 

fst is the direct tension strength; 

a, h are the notch depth and the beam height, 
respectively; and 

Gr is the fracture energy. 


3.2 Finite element method (FEM) 


The computational analysis was carried out with software 
ANSYS [6] and choosing its element SOLID185 to model 
the concrete in 3D. 


After the concrete experiences a cracking phase, the 
internal forces are transmitted to the fibers, which then 
govern the behavior of the material. 


The Multilinear Material Model used in this work 
(CAST) can approximate behavior laws both in 
compression and in tension. Fig. 4 describes the boundary 
conditions of the beam considered in the model. 
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The UHPFRC was simulated as a composed material 
with a law in compression that was obtained from 
experimental data and a law in tension from an Inverse 
Analysis that includes the material’s post-cracking 
behavior. 


SOLID185 is a 3D element that allows considerations 
to represent plasticity, hyperelasticity, large displacements, 
and large strains. It also allows simulations of quasi- 
incompressible elastoplastic materials and fully 


incompressible hyperelastic materials. 


The element is defined by eight nodes with three 
degrees of freedom each (translations in x, y, and z 
directions), as shown in Fig. 4. 


Fig. 4: Boundary element and element SOLID 185. 


CAST is an elastic isotropic multilinear material with 
the same elastic behavior in compression and in tension, 
but with elastic limit and isotropic hardening behavior that 
can be different in each case. The behavior in tension uses 
the Rankine criterion, while the behavior in compression 
uses Von Mises. 


The UHPFRC properties, such as its modulus of 
elasticity and its Poisson’s coefficient, had to be known for 
the simulations. These values were maintained constant in 
each specimen that was modeled. 


The behavior laws in tension and in compression were 
different in each specimen since those behaviors were 
drawn from the experimental tests and the results from the 
Inverse Analysis. 


3.3 Forces in cracked section 


Fig. 5 shows the cracked cross-section of a prismatic beam 
subjected to bending forces, and where two different 
regions can be easily identified. 


Firstly, there is the zone without any cracking, which is 
the part of the section where the force distribution 
corresponds to a linear elastic behavior. 


Secondly, there is the cracked zone, which is the part 
of the section where the force distribution directly depends 
on the effectively of the fibers inside the concrete matrix, 
which can be determined via Inverse Analysis. 
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Fig. 5: Forces in the cracked section. AFGC [4] modified. 


The force equilibrium in the section results in 
equations 4 to 8 with “b” identifying the contribution of 
the regions with cracks, while “f” identifies the cracked 
ones. 


Nbe= = .ob.c.b (4) 
Nbt == .ot.c’.b (5) 
Nb = Nbc + Nbt (6) 
M=Mb+Mf (7) 
N=Nb+Nf (8) 


The system of eight equations to solve is bound to 
equations 9 to 16, shown in the following. 


1 
Nb == -E.Xm.b.h?.[4?—B?] (9) 
Nf, =N}. C. D+K.b.H.(1—D) (10) 
1 
Mb=3 .E.Xm.b.h*.[A?— B?]+h.an.Nb (11) 


Mf, = Mfj.(C.D)? + K N$£,4.0 —L.O7of,,, (12) 


N = Next = Nb+Nf (13) 
M = Mext = Mb + Mf (14) 
ot=E.Xm. h.(an — a) (15) 


= |x 2 16 
w= m+2 g| 3 (16) 


where: 

h.a is the relative length of the crack; 

h.an is the relative height of the neutral axis; 
Xm is the curvature of the region without cracks; 


b, h are the width and height of the beam cross-section, 
respectively; 


I is the moment of inertia of the rectangular section; 
and the variables are: 
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A=1-«,; B=x—«,; C= p= O 
ij Gisa 
A A (ipa h)* b 
H = of; ae K — Cipt h; L = Sir ) : 
0=1i-D 


The equation system is solved for each point of the (o 
vs. œ) curve by using the known experimental points (F, 
@) and the parameters calculated in the previous step. 


3.4 Validation using energy fracture 


The area under the analytical (o vs. œ) curve, which is 
obtained via the Inverse Analysis commented in the 
previous section, represents the fracture energy, GF, of the 
material. In the same form, the area under the experimental 
(F vs. ò) curve gives a measure of GF, calculated 
according to specifications given by RILEM TC50 [5]. 


The Fracture Energy can be found using the load- 
displacement data and the equation 17. A graphical 
comparison is made between both behaviors and the 
fracture energy is then calculated for every specimen with 
1% and 2% of fiber content. 

H; 


Gr = I-a (17) 


where: 


Wr is the total area of the curve under the graphic of 
load versus deflection; 


b is the thickness of the beam (mm); 
h is the height (mm); and 


a is the length of the notch made in the lower center of 
the beam. 


3.5 Validation using finite element method 


The law of behavior in compression is obtained from the 
experimental data, and the law of behavior in tensile is 
obtained by inverse analysis. 


The numerical simulation of the flexural test is 
performed using these behavior curves as input data. An 
analytical load-displacement curve is obtained for each of 
the specimens. Then this analytical curve is compared with 
the response obtained experimentally, as initially indicated 
in Fig. 1. 

The approximation between the analytical and 
experimental curves, indicated above, is a measure 
adopted in this investigation to validate the inverse 
analysis. With this, it is possible to verify the effectiveness 
of the methodology proposed by the French standard in the 
AFGC [4], developed from the mechanical equilibrium of 
Fig. 5 and by equations 4 to 16. 
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IV. EXPERIMENTAL AND NUMERICAL RESULTS 


4.1 Compressive strength and modulus of elasticity 


Twenty cylinders measuring 5cm (2in.) in diameter and 
10cm (4in.) in length were used for uniaxial compressive 
tests for the material with 1% of fibers. The compressive 
strength was calculated as the average for those twenty 
specimens, resulting in 151MPa (22ksi), as shown in 
TABLE 3 (fcm=151MPa). The standard deviation was 
4.3MPa (0.7ksi). 


Table 3: UHPFRC compressive strength and modulus of 
elasticity (IMPa=145psi). 


Specimen œ (MPa) E (MPa) 
1 154 50709 
2 147 44504 
3 146 46104 
4 158 49551 
5 150 48209 
6 150 46780 
7 146 47503 
8 152 45802 
9 153 44548 
10 155 46556 
11 144 43768 
12 159 50799 
13 147 47035 
14 150 49184 
15 150 49129 
16 150 48193 
17 146 45522 
18 150 49293 
19 152 49595 
20 158 51374 

Average: 151 47708 


In each test, the o-s curve is obtained and the modulus 
of elasticity is calculated by a linear approximation 
between 5% and 80% of the compressive strength, which 
averaged 47708MPa (6919ksi), as also shown in TABLE 
3. The standard deviation was 2.2MPa (0.3ksi). 


Fig. 6 shows behavior the average curve in uniaxial 
compression for the tested specimens. The maximum 
compressive stress (fem) value occurs for a strain value 
equal to 0.0033. The characteristic resistance value (fck) 
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was 143.43MPa (20.8ksi) with a 95% probability of 
exceedance, obtained using the Student-fisher law. 


Cylinders containing 2% of fibers were tested; in all 
cases, the results showed resistance values significantly 
lower than those manufactured with 1% of fibers and 
therefore were discarded. We presumed that reduction of 
strength was due to fiber agglomerations and the formation 
of internal voids. 


c (MPa) 


— Average experimental curve 


0 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 


z (mm/mm) 


Fig. 6: UHPFRC compression constitutive behaviour. 


The description in the previous paragraph meets the 
AFGC [4] recommendations, which proposes to 
characterize the compression behavior of UHPFRC 
according to the values of the characteristic compression 
strength and the modulus of elasticity. 


4.2 Constitutive behaviour in tensile 


Fig. 7 and Fig. 8 show (F vs. 6) curves for each of the 
tested beams, as well as the average curve and considering 
fiber content of 1% and 2%, respectively. 


— Average experimental curve (1% fiber) 


Fig. 7: Experimental curves for beams with 1% of fibers in 
bending (IkN=225lbf; 1mm=0,04in). 


The area under each curve was calculated to determine 
the fracture energy according to RILEM TCS0 [5]. 
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The values of the elastic load, Fte; of the elastic 
strength in tension, o; and of the deflection, dte; are 
presented in TABLE 4. 


Average experimental curve (2% fiber) 
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Fig. 8: Experimental curves for beams with 2% of fibers in 
bending (1kN=225lbf; Imm=0,04in). 


Table 4: UHPFRC elastic load and deflection, and 
strength in bending (IkKN=225lbf; Imm=0.04in; 
IMPa=145psi) 


Fte (kN) dte (mm) o (MPa) 


Specimen 1% 2% 1% 2% 1% 2% 


CP-1 9.9 20.0 0.032 0.030 9.1 18.4 
CP-2 10.0 15.1 0.039 0.024 9.2 18.5 
CP-3 10.1 15.1 0.039 0.030 9.3 13.8 
CP-4 10.6 10.3 0.037 0.020 9.7 9.5 
CP-5 10.2 0.018 9.3 
CP-6 10.3 0.021 9.4 


Average 10.1 13.5 0.037 0.024 9.3 13.2 


TABLE 5 presents the results obtained from the post- 
cracking behavior for the maximum load Ftcr and its 
corresponding deflection ôtcr. 


Table 5: UHPFRC inelastic load and deflection in bending 
(1kN=2251bf; Imm=0.04in) 


Ftcr (kN) 

Specimer 1% 2% 1% 2% 
CP-1 9.2 21.0 0.83 0.81 
CP-2 11.4 20.8 1.30 1.03 
CP-3 13.8 22.5 1.04 1.08 
CP-4 14.3 22.0 1.28 1.04 
CP-5 23.4 0.91 
CP-6 26.4 1.09 
Average 12.2 22.7 1.11 0.99 


dtcr (mm) 
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Fig. 9 and Fig. 10 show (o vs. @) curves obtained from 
Inverse Analysis for each one of the beams, with fiber 
contents of 1% and 2%, respectively. 


The area under each curve was calculated to determine 
the fracture energy. 
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Fig. 9: Numerical (o vs. œ) curves for beams with 1% of 
fibers in bending (IMPa = 145psi; Imm = 0.04in). 
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Fig. 10: Numerical (o vs. œ) curves for beams with 2% of 
fibers in bending (IMPa = 145psi; Imm = 0.04in). 


It was calculated from the relations (F vs. 6) and (o vs. 
@) as is showed in TABLE 6 for each one of the tested 
specimens, where a good fit can be observed between the 
two averaged results. 


Table 6: Fracture Energy (GF) for UHPFRC beams with 
1% and 2% of fibers (1kJ/m2 = 0.0006BTU/in2). 


Fracture Energy (kJ/m?) 


1% of fibers 2% of fibers 
Inverse RILEM Inverse RILEM 
Analysis TC50- Analysis TC50- 
AFGC FMC AFGC FMC 
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CP-1 16.91 11.66 26.53 23.76 
CP-2 23.71 18.75 29.29 19.17 
CP-3 23.70 15.42 29.19 23.10 
CP-4 23.32 19.22 18.87 21.76 
CP-5 3249 32.40 
CP-6 32.48 57.96 
Average 21.91 16.26 26.67 24.04 


Fig. 11 and Fig. 12 show (o vs. £) curves obtained from 
the transformation of q into e using equations 1 to 3, with 
fiber contents of 1% and 2%, respectively. 
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Fig. 11: Numerical o vs. e curves for beams with 1% of 
fibers in bending (IMPa = 145psi; Imm/mm = lin/in). 
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Fig. 12: Numerical o vs. e curves for beams with 2% of 
fibers in bending (IMPa = 145psi; Imm/mm = Tin/in). 


Fig. 13 and Fig. 14 show the results of computational 
modeling, with fiber contents of 1% and 2%, respectively. 
Models and experiments showed good agreement. 
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The behavior analytic curves obtained in this research 
showed a similar trend with the curves (F vs. 5) obtained 
by Denairé et al. [14] using inverse analysis. Similarly, 
Mezquida et al. [15] carried out inverse analysis 
methodologies, based on the closed-form non-linear hinge 
model, to define the material's behavior. They obtained a 
similar response to this research in both cases: when the 
UHPFRC exhibits strain-hardening constitutive stress- 
strain behavior and when it exhibits strain-softening 
behavior. 


apo ces Experimental response with 1% fibers 
144 Analytical response with 1% fibers 
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Fig. 13: Average response for 1% of fibers. Experimental 
vs. numerical (IkN = 2251bf; Imm = 0.04in). 


Arenes Experimental response with 2% fibers 
284 Analytical response with 2% fibers 
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Fig. 14: Average response for 2% of fibers. Experimental 
vs. numerical (IkN = 2251bf; Imm = 0.04in). 


Also, Chanvillard and Rigaud [16] studied three points 
bend test on notched specimens and applied an inverse 
analysis to extract the tensile strength versus crack opening 
relationship. Again, the behavior curves showed a similar 
trend to the results in this research. 
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V. CONCLUSIONS 


The results obtained from this investigation allow the 
following conclusions: 


The computational modeling of UHPFRC beams can 
be satisfactorily carried out by considering the behavior of 
the composite material under a homogeneous premise. 
This can be accomplished with bending tests and the 
determination of behavior laws for the matrix with fibers 
in uniaxial compression and tension; 


The constitutive laws for the UHPFRC material were 
experimentally and numerically determined for each of the 
beams considered. The o vs. € curves obtained in each case 
were considered as input data for the computational 
modeling carried out in Finite Elements in ANSYS. The 
results generated numerical F vs. 6 curves that were 
compared with the ones experimentally obtained, showing 
a good fit between them; 


The finite element SOLID185 used to model the 
matrix, together with the CAST material model used to 
simulate the behavior of the cracked section governed by 
fibers, were adequate to model the UHPFRC; 


The Inverse Analysis procedure showed to be adequate 
to determine the behavior curve in tension of the 
considered beams made of UHPFRC, even considering the 
post-cracking response of the material; 


The validation of the Inverse Analysis by means of 
calculating the fracture energy showed to be satisfactory 
for the beams with 2% of fiber content. The average value 
calculated from (o vs. œ) numerical curves was 27 kJ/m2 
(0.017BTU/in2), while the value obtained from the 
experimental (F vs. 6) curves was 24 kJ/m2 
(0.015BTU/in2), i.e., a difference of roughly 10%. 
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